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Abstract 

That muon neutrinos oscillating into the mixture of tau neutrinos I'r and sterile neutrinos i^s 
has been studied to explain the atmospheric disappearance. In this scenario, the effect of Earth 
matter is a key to determine the fraction of Vg- Considering that the Earth matter density has 
uncertainty and this uncertainty has significant effects in some neutrino oscillation cases, such as 
the CP violation in very long baseline neutrino oscillations and the day-night asymmetry for solar 
neutrinos, we study the effects caused by this uncertainty in the above atmospheric oscillation 
scenario. We find that this uncertainty seems to have no significant effects and that the previous 
fitting results need not to be modified fortunately. 
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To explain the atmospheric muon neutrinos z/^ disappearance, the scenario of z/^ oscillat- 
ing into Z/+ has been studied 0,0], where z/+ is the mixture of tau neutrinos and sterile 
neutrinos z/^, and defined as Vj^ = VrCos^ + VgSin^. Since for so-called "matter effects" |4|, 
the oscillation probabilities P^^^^^ and P^^^^^ are different for a muon neutrino with certain 
energy that travels a distance in Earth, one can expect togive a limit on ^. The reported 



ma. 



results from Super-kamiokande have given limits on sm^^[ 

In the calculation of Pv^^^vs ■, the neutron number density of Earth is a critical quantity. 
However, today the knowledge of Earth matter density which determines the neutron num- 
ber density is only to some certain precision ^. As to the preliminary reference Earth model 
(PREM) 



TV IS only to some certam precision pj. As to tne prenmmary reierence Hia:. 
a], the uncertainties due to the local variation have been documented^]- Quanti- 



tatively its precision is roughly 5% averaged per spherical shell with thickness of 100 km or 
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The effects of Earth matter density uncertainty have been studied in some neutrino 
oscillation cases, such as the CP violation in very long baseline neutrino oscillations 0, 
llo| and the day-night asymmetry for solar neutrinos [11]. Ones find this uncertainty has 
significant effects in these casesp, Q, U]. Since the Earth matter is a key to determine 
the fraction of z/^, this uncertainty could also have an effect on the limit of sin^^. In this 
brief report, we study the density uncertainty in Earth matter and then investigate its 
implications on the results of siri^S,- 

We begin our discussions with the effective Hamiltonian that governs the propagation of 
the neutrinos in matter. In the (2+2) mo dels jl^. the relevant (z/^, z/_|_) evolution is given by 
the Schrodinger equation 





~ H(x) 











(1) 



with the effective Hamiltonian I 



, , Am2 / -cos2^ sin2^ \ ^ / \ 

H{x) = —— + V2Gf . (2) 

y sin 2^ cos 2^ J y i sin^ ^A/'„(x) J 

Here is the neutrino energy, Am^ and 6 are the usual mass and mixing parameters in the 
z/^ Ut- oscillating model, Gp is the Fermi constant, and Nn{x) = p{x)Na{^ — ydx)) is the 
neutron number density with p{x) the matter density in g/cm^, A^^ the Avogadro number 
and Ye{x) the electron number fraction, respectively. For increasing values of sin^E,-, we get 



a smooth interpolation from i/^ oscillations {si'n?^ = 0) to pure z/^ Ug oscillations 

{sin'^C, = 1), passing through mixed active-sterile transitions (0 < sin'^C, < !)• Replacing 
Nn{x) by —Nn{x), we can also get the effective Hamiltonian for relevant antineutrinos. 

Now we consider the uncertainty in Earth matter and its implications on the atmospheric 
neutrino oscillations. In some simple cases, for example, if the neutron number density Nn{x) 
suffers from a global shift (independent of x), the induced effects on sin^^ and Nn{x) are 
degenerate in the effective Hamiltonian Q , such as sin^ (i±5%) ^^^^ ^ when Nn 
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{l±5%)Nn{x). 

Generally, at a given point x in Earth, the available matter density, which determines 
the neutron number density, is an average value with some prescribed errors, such as the 
widely used PREM model 6|. We can define the average density p{x) as an average over all 
samples of density profiles {p{x)} 

p{x) ^ {p{x)) = j V[p{x)]F[p{x)]p{x) , (3) 

and the error cr(x) as a variance function 

a{x) ^ ^{p^x)) - , (4) 

where F[p(x)] is the probability density of the density sample p{x). Accordingly, the av- 
eraged probability for the a flavor neutrino oscillating into the j3 flavor neutrino should 
be 

(P.,^.,(L,E,)) ^ j V[p{x)]F[p{x)]P,^^.^mE,,p{x),Y,{x)) . (5) 

with L the neutrino's travelling distance in Earth. Furthermore, we can write the variance 

as 

<5P,„^,,(L,E,) ^^j V[p{x)]F[p{x)]{P,^^.^{L,E,,p{x),Y,{x)) - (P,„^,,(L, E,)))^ . (6) 

In this brief report, we introduce a logarithmic normal distribution^^ to represent the 
probability density function of the Earth matter density samples 

s{x) = y'ln[l + r2(x)] , (8) 
Pq{x) = p{x) exp[-^s2(x)] , (9) 



where r{x) = (j{x)/p{x) characterizes the precision of Earth matter density. And then we 
use Monte Carlo calculations to generate the values of F[p{x)] between and 1 at a given 
point X along the propagating path of the neutrinos. With the chosen p{x) and r{x), we 
obtain the value of p{x) from Eq.((7j) by computing the values of s{x) in Eqs.® and po{x) in 
Eq.Q. Hence, the averaged oscillation probability (0) and the corresponding variance (jH)) 
can be calculated. Specifically, we take Ye{x) = 0.5 and r(x) = 5% with p(x) given by the 
PREM in our numerical calculations. 

Usually, the experiment results are reported as the event number, which can be calculated 

as 

N = N,ct>,SEu){Pu^^u,m E,))ap{E,)CdE, 
Je, 



= E No<P.^{Et){P.^^.,{L, Et))a^{Et)CAEt . (10) 
fc=i 

where E^, Ei denote the upper bound and lower limit of the detecting energy, Nq is a 
normalization factor with unit conversions, (pu^ is the a flavor neutrino beam flux spectrum, 
is the charged current cross section of /? flavor neutrino, C is the product of the detector's 
size and running time, and AE^^ is the kth energy bin size. Accordingly, we can define 



V 



6N = / No(j),AEu)SP,^^,,{L, E,)a^{E,)CdE, 

J El 



= Y: Noct^.SEt)5P.^^.,{^ Et)ap{El)CAEt , (11) 
fc=i 

as the variance and relative variance of event number caused by the uncertainty in Earth 
matter density, respectively. 

For example, in the numerical calculations we take the relevant data listed in the table 

Ti n 

II, IV, V of [14] and table III of [15|]. We use the mass and mixing parameters: Am = 



2.0x lO-^eV^ sin^ 26 = 1.0, sin^ ^ = 0.19|3|, and neutrino travelling distance L = /-2r cqs_e, 
where / is the neutrino production height in atmosphere (slant distance in km, listed in 1^), 
r = 6371km is the Earth radius and cosG is the zenith angle. We plot the event number 
produced by the atmospheric muon neutrinos and antineutrinos in six zenith angles with 
three energy bins as well as their variances arising from the uncertainty of Earth matter 
density by using Eq.()10|) and (lllj). As shown in Fig.l, we find the variance of event number 
is smaller than ~ 10^^, accordingly the relative variance is smaller than ~ 10~^, for the 



longest baseline cos = —1.0 with the largest matter effect, while in the other baselines, the 
variance and also the relative variance are much smaller. Hence we can draw a conclusion 
that the uncertainty of Earth matter density seems to have no significant effects in the 
oscillating scenario of z/^ — > z/^cos^ + UgSinC,- 

In summary, considering that the Earth matter density has uncertainty and this un- 
certainty has significant effects in some neutrino oscillation cases, such as the CP vio- 
lation in very long baseline neutrino oscillations and the day-night asymmetry for solar 
neutrinos, we study this uncertainty in the atmospheric neutrino oscillating scenario of 
J^fi — ^ Ut-cos^ + UsSin^, and analyze the effects caused by this uncertainty on the previous 
fitting results. We find that this uncertainty seems to have no significant effects and need 
not to modify the previous fitting results fortunately. 
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FIG. 1: Event numbers plotted as a distribution with zenith angle cos0 in three neutrino energy 
bins. The square is Honda^J] expectation for no oscillations. The dot is given by substituting 
{P) for P in Eq.(8) and up-triangle represents the variance defined in Eq.(9). These plots use the 
mass and mixing parameters: Am^ = 2.0 x lO^^eV^, sin^ 20 = 1.0, sin^ ^ = O.IQ^], and neutrino 
travelling distance L = 1 — 2rcos0, where / is the neutrino production height in atmosphere(slant 
distance in km, listed in^^) and r = 6371km is the Earth radius. 
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